function sys=plant(t,x,U1,U2,D)
sys=zeros(4,1);

sys(1)=x(2);  
sys(2)=U1+D;
sys(3)=x(4);  
sys(4)=U2+D;
end